


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations l. Thesis and Dissertation Collection, all items 


1977-09 


Experimental study of natural convection in 
thin horizontal liquid layers 


Yurteri, Mahmut 


http://hdl.handle.net/10945/18123 
Copyright is reserved by the copyright owner. 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
D U DLEY research materials and institutional publications created by the NPS community. 
qe Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
KNOX appointed — and published — scholarly author. 


LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http: //wwwenps.edu/library 


EXPERIMENTAL STUDY OF NATURAL CONVECTION IN 
THIN HORIZONTAL LIQUID LAYERS 


Mahmut Yurteri 








NAVAL POSTGRADUATE SCHOOL 


Monterey, California 





TRESS 


EXPERIMENTAL STUDY OF NATURAL CONVECTION IN 
ПЕ ТООЛОН LIQUID LAYERS 


Бу 


Mahmut Yurteri 


September 1977 


Thesis Advisor: M.D. Kelleher 





Approved ror public release; distribution unlimited. 


1182126 





UNCLASSIFIED 


SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 


REPORT DOCUMENTATION PAGE 


1. REPORT NUMBER 2. GOVT ACCESSION NO. 


4. TITLE (and Subtitie) ۱ 
Experimental Study of Natural Convection 


m тп Herizontal Lrq@uid Layers 








READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


5: ТҮРЕ OF REPORT & PERIOO COVEREO 
Engineer's Thesis; 
September 1977 


6. PERFORMING ORG. REPORT NUMBER 
















7. AUTHOR(e) 8. CONTRACT OR GRANT NUMBERY(4) 


Mahmut Yurteri 


9. PERFORMING ORGANIZATION NAME AND ADORESS 10. PROGRAM ELEMENT, PROJECT, TASK 


AREA & WORK UNIT NUMBERS 
Naval Postgraduate School 
Monterey, California 93940 


11. CONTROLLING OFFICE NAME AND AOORESS 12. REPORT OATE 
Naval Postgraduate School September 1977 


Monterey, California 93940 13. m OF PAGES 


4. MONITORING AGENCY NAME & AOORESS(i(f different (rom Controlling Office) 1$. SECURITY CLASS. (ot thia report) 


Unclassified 


15а. OECLASSIFICATION/ DOWNGRADING 
SCHEDULE 





16. DISTRIBUTION STATEMEN T (ol thie Report) 


Approved for public release; distribution unlimited. 


17. DISTRIBUTION STATEMENT (ol the abatract entered in Block 20, if dií(erent from Report) 


18. SUPPL EMENTARY NOTES 


19. KEY WORDS (Continue on reverse side if necessary and identify by block number) 


Natural Convection 

Heat Transfer 

Thin Liquid Layers 
Critical Rayleigh Number 


20. ABSTRACT (Continue on reverse aide if neceeeary and identify by biock number) 


Heat transfer through a very thin horizontal liquid layer 
bounded on top by a cooled glass plate and on bottom by a 
heated copper plate was measured. Results were correlated 
in terms of Nusselt number as a function of Rayleigh number. 
Three different liquids were used under conditions to give 
a range of Rayleigh number from 350 to 4100 and Prandtl number 





FORM 
рр , 5295, 1473 EDITION OF 1 NOV 65 15 OBSOLETE 
(Page 1) S/N 0102-014- 6601 | _ _ UNCLASSIFIED  — . 

l SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 





UNCLASSIFIED 


$.cuwT Y CLASSIFICATION OF THIS PAGESWHen Dora Enteros 





(20. ABSTRACT Continued) 


From 34 to 477. 
Up to a critical condition heat is transferred by 


conduction alone.  Convection appears when Ra = 1600 + 100. 
A correlation for critical Rayleigh number as a function of 


Prandtl number was obtained. 


1473 UNCLASSIFIED 


I 
Жо, l nie 
3 N 0102-014-6601 2 SECURITY CLASSIFICATION OF THIS PAGE(When Deta Entered) 





Approved for public release; distribution unlimited. 


Experimental Study of Natural Convection in 
Thin Horizontal Liquid Layers 


by 


Mahmut Yurteri 
Lieutenant, furkish Navy 


Submitted in partial fulfillment of the 
requirements for the degrees of 


MASTER OF SCIENCE IN MECHANICAL ENGINEERING 
and 


MECHANICAL ENGINEER 


кош the 


NAVAL POSTGRADUATE SCHOOL 
September 1977 





ABSTRACT 


Heat transfer through a very thin horizontal liquid 
layer bounded on top by a cooled glass plate and on bottom 
by a heated copper plate was measured. Results were 
correlated in terms of Nusselt number as a function of 
Rayleigh number. Three different liauids were used under 


conditions to give a range of Rayleigh number from 350 to 


٦٢ and Prandtl number Irom 34 to 477.‏ س۲۰۰۰ 


Up to a critical condition heat is transferred by 
conduction alone. Convection appears when Ra = 1600 + 100. 
A correlation for critical Rayleigh number as a function of 


Prandtl number was obtained. 
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А. BACKGROUND 

The field of natural convection in enclosures encom- 
passes many variations of geometry and input situations. 

One of these situations is heat flow in the narrow space 
between two horizontal parallel plates where the lower 

plate has a higher temperature than the upper plate. When 
the lower plate is heated, the fluid remains immobile and 

an unstable stratification is formed, inasmuch as the warmer 
fluid of lower density is located below the cooler fluid 
whose density is higher. 

In his classical work on convection currents ina 
horizontal layer of fluid Rayleigh [Ref. 1] examined the 
case of temperature gradients in a layer of fluid. He 
tormulated the flow equations for a discrete disturbance 
in the fluid and determined the conditions under which the 
disturbance would amplify causing the layer to become 
unstable. Rayleigh recognized that the unstable stratifi- 
cation must break down at a certain value of the temperature 
difference above which a convective motion must be generated. 


Below the critical value of Rayleigh number, defined as: 


5 
E n 


a V 


pure conduction is observed. Above the critical value of 


the Rayleigh number convection begins. 
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Benard [Ref. 2] performed the first laboratory controlled 
experiments on thermally unstable liquid layers. He worked 
with liquid layers on the order of several millimeters, 
lying on a metallic plate which was heated and maintained 
at a uniform temperature. The upper surface of the liquid 
layer was free and at a lower temperature than at the plate 
surface, since it was in contact with the ambient air. 

Benard described the change from conduction to convection 

in two phases. The fluid remains at rest until the vertical 
temperature difference becomes sufficiently large. A pre- 
liminary motion of the fluid then results. Shortly thereafter, 
this first phase of relatively short duration appears, in 
which the fluid forms cells of nearly regular polygons with 
four to seven sides. This phase lasts from a few seconds 

to several minutes for more viscous fluids. During the 

second phase, the cells become equal and reqular and align 
themselves. The limit of the second phase is thus a permanent 
regime of cells with vertical boundaries and hexagonal 

Eros s sections as’shown in Figure I. Liquid rises in the 

core of the cell, moves outward at the top, descends at 

the vertical boundary between adjacent cells and moves 

inward at the bottom. 

In 1926 Jeffrevs (Ref. 3] presented a solution to 
Rayleigh's problem for two boundaries with no tangential 
forces between boundary and fluid. His study included 


solutions for the case of two rigid conducting boundaries 


5 


- 


Е 
u 

s 

= E 7 mE 
u u 


SSW ath, 


873 тї час 





O 





О 


oes o o 


View 


Plan 





Section 


Vertical 


Sehem Skereh of Benard Cell Structure 


500058 l. 





and a rigid conducting boundary at the base with a free 
surface at the top. For two rigid conducting boundaries 
he found the critical Rayleigh number to be about 1200. 
By a revised solution method in 1928 [Ref. 4] Jeffreys 
obtained 1709.5 for the critical Rayleigh number for two 
rigid conducting boundaries. 

ег Low [Re£. 5], Sutton [Ref. 6], Malkus [Ref. 7] 
and Catton [Ref. 8] in their theoretical works agreed on 
Ee critical Rayleigh number of 1705 + 5. 

Schmidt and Milverton [Ref. 9] with their admittedly 
rough experiments with distilled water determined the 
critical value of the Rayleigh number to be 1770. In 1938 
Chandra [Ref. 10] observed a columnar motion well under the 
critical Rayleigh number. His test fluid was alr with fluid 
layer thicknesses varying between 4 and 16 millimeters. 

He concluded that for all fluid layer thicknesses below 

10 mm. a type of motion other than cellular motion which 
he called columnar motion occurs below the critical Rayleigh 
number. Later in 1952 DeGraaf and Van Der Held [Ref. 11] 
worked with horizontal and inclined air layers between 
parallel plates. They found that in horizontal air layers 
the air remains at rest until the Rayleigh number reaches 
about 2000. When this value is exceeded, the cellular 
motion sets in, first in the form of hexagonal prisms, but 
with a tendency to change into rows of tetragonal prisms 
when the Rayleigh number increases. In contrast with the 


Benard cells in liquids, the air descends in the middle. 
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This was explained by the increase of the kinematic 
viscosity with increase in temperature for gases. At 
smaller fluid layer thicknesses of 5.5 and 6.9 mm. a 
columnar motion, first described by Chandra, was observed 
when the Rayleigh number exceeded 1400. For these thick- 
nesses, cellular motion was observed for Rayleigh numbers 
above 1600. This low critical Rayleigh number for thinner 
air layers was explained by the authors as being due to the 
very great temperature gradients. 

Sutton [Ref. 6] explained previous works in his 1950 
paper. He showed that the criterion for the "columnar" 
mode involves only the ratio of the absolute temperatures 
of the upper and lower surfaces. He also derived an 
expression and verified the critical temperature difference 
at which the transition takes place from the "columnar" 
to "cellular" mode. He also showed that the "cellular" 
mode will occur if the depth of the test fluid exceeds a 
certain value, but that for the more shallow layers, the 
"columnar" mode will be generated initially, ultimately 
passing to the "cellular" mode for increased temperature 
difference. 

Malkus [Ref. 7] performed a series of experiments with 
distilled water and acetone. The distance between the two 
parallel plates ranged from 0.05 inches to 8.0 inches during 
these experiments. His experimental critical Rayleigh 


number was found to be 1700 + 80. 


l6 


( 





Ernst Schmidt and Silveston [Ref. 12] examined heat 
transfer through a horizontal liquid layer bounded on the 
top by a cold surface and on bottom by a heated surface. 
Five different liquids were used at different temperatures 
for a range in Prandtl numbers of from 3 to 4000. Test 
liquids were distilled water, ethylene glycol, heptane and 
two silicone oils. Layer depths ranged from 1.45 mm. to 13 
mm. Optical observations were made on the patterns formed 
in the layer by convection. They concluded that up to a 
critical condition heat is transferred by conduction alone 


and convection appears at 


каға 170005 3%; 


They observed three distinct convection regimes. The first 
occurred as convection sets in and appeared to be a honey- 
comb pattern in the layer. As the Rayleigh number increases, 
the pattern was found to change to a series of stripes 

where the heat transfer had laminar character. The third 
regime was observed at higher values of Rayleigh number 

with a tangled and disordered pattern which was the turbulent 
region. Results, including the data of Mull and Reiher 

[Ref. 13] were correlated by Nu as a function of Ra. The 
Manors comeluded that the data indicated four distinct 
modes of heat transfer. Above Ra = 10,000 they found that 
Rayleigh number is not the only correlating parameter, and 


the Prandtl number must be taken into account. 
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The correlations derived from the data were: 


Creeping region: Ып = 0 0012 (Ra) 9:90 
Laminar region: Nu = 0.24 (Ra) 2> 
Transition region: Nu = 0.30(Gr) 9: 19 (p, 0: 21 
= 0 ж (Ба)? ° (ру) 06 
КОО ШО region: Nu= 0,10(6r) 1 (ре) 929 
ری ۱۰۰ ۷ا2‎ ^+ 


The authors commented that the creeping convection 
begins about Ra = 1700 and the laminar region about 
Ra = 3000. Starting point of the transition region was to 
be determined by Ra = 8000 (ру) 9°“ whereas the starting point 
of the turbulent region was given by Ra - 18,000 Pr 

Heat transfer measurements in horizontal fluid layers 
seem to be generally lacking in the literature except for 
some extensive data on air obtained by Mull and Reiher 
(Ref. 13] and on liquids by Ernst Schmidt and Silveston 
[Ref. 12]. One of the objectives of this study was to 


obtain data for liquids near the critical Rayleigh number, 


which is the area least covered by the previous works. 


ES THESIS OBJECTIVES 
The objectives of this study were to investigate natural 
convection in thin horizontal liquid layers heated from 


below, in particular to investigate the change from the 
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conduction heat transfer regime to the convection heat 
transfer regime. It was intended to determine the critical 
Rayleigh number associated with the transition. 

Another objective of the study was to determine corre- 
lations of the Nusselt number as a function of the Grashof 
and the Prandtl numbers. 

To accomplish these objectives a series of experiments 
with several different fluids were conducted. An experimental 
apparatus was designed where the hypothetical system of a 
system fluid contained between two infinite, horizontal 
surfaces was closely approximated by containing the fluid 
between two parallel plates and using layer depths which 


were very small compared with the dimensions of the plates. 


19 





II. APPARATUS 


A. DESIGN CONSIDERATIONS 

The hypothetical system of a fluid contained between 
two infinite, horizontal, conducting surfaces was closely 
approximated by containing the test liquid between parallel 
Square plates and using layer depths which were very small 
compared with the dimensions of the plates. Larger layer 
depths were avoided to keep the edge effects to a minimum. 

Several preliminary experiments were performed with 
distilled water to obtain an appreciation for the general 
performance of the apparatus. These preliminary experiments 
gave a general idea of the required properties of the test 
liquid and the required distance between the two plates in 
order to obtain Rayleigh numbers about the critical value 
of the Rayleigh number. The data obtained from these pre- 
liminary experiments are not included in this study because 
the obtained Rayleigh numbers were outside the range of 
interest. 

A petroleum based oil, commercially known as Mobil 603 
was chosen as the second test liquid after these preliminary 
experiments. Two sets of data were obtained with Mobil 
603. A plain glass plate was used as the top surface of the 
test chamber during the first set of experiments. During 
these experiments a problem with short circuiting of the 


heater with the copper plate was encountered at high 
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temperatures. After these experiments a very thin, heat 
resistant gasket was put between the copper plate and the 
electric heater to increase the electrical insulation. 
During the second set of experiments a glass plate coated 
with a transparent, electrically conductive coating was 
used. The data obtained from these two sets of experiments 
are contained in the study. The data of the second set of 
experiments was consistent with the data obtained later and 
this data is used for correlations. 

With the help of the preliminarv water data a charac- 
teristic plate spacing L, was selected to obtain a Rayleigh 
number close to the critical Rayleigh number. Rayleigh 
numbers above and below the critical value were obtained 
by changing other parameters such as the temperature 
difference between the two plates of the test chamber and 


the fluid properties. 


B. DESCRIPTION OF THE APPARATUS 

The test apparatus as shown schematically in Figure 2, 
consisted of three principal components. These were a test 
chamber where the test liquid was contained between two 
parallel plates, a cooling chamber above it where cooling 
water circulated and a quard heater/insulation assembly. 

A copper plate with an attached heater at the bottom 
and a glass plate at the top with a plexiglass shim between 
constituted the test chamber. The cooling chamber consisted 


of a plexiglass plate on top, the glass plate of the test 
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chamber on bottom and a thick plexiadlass spacer with inlet 
and outlet tubes on the sides. Layers of insulation with 
six thermocouples in two groups and a guard heater between 
the insulation layers constituted the insulation piece. 
The guard heater and the thermocouples were placed in the 
insulation to minimize the heat loss. 

A one-quarter inch thick and seven inch square copper 
olate, as shown in Figure 3, was used as the lower heated 
surface of the test chamber. Both sides were given a 
machine finish and later the upper surface was polished to 
a mirror finish. Flattened beads of five thermocouples 
(30 gage copper-constantan) were attached on this surface. 
The thermocouples were painted with a nonconducting paint 
E vosdoshortcircuiting through the copper plate. The 
thermocouple leads were taken out through 0.05 inch diameter 
holes bored into the copper plate. These two holes were 
sealed with silicon rubber. 

A third channel was bored into the copper plate and a 
small diameter stainless steel tube was soldered to the 
outer end of it. A plastic tube connected this metal tube 
to a container filled with test liquid which acted as an 
expansion tank for the test chamber. It also kept the test 
chamber near atmospheric pressure. 

The heater element was made of twelve feet of nichrome 
wire having a total resistance of twelve ohms. The twelve 
foot long wire was shaped into parallel strips, one-quarter 


inch apart covering an area six inches square as shown in 
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FIGURE 3. Location of the Thermocouples on the 
08) 
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. Figure 4. The shape of the heater element and the thickness 
of the copper plate insured a uniform heat flux. The heater 
element was attached to the bottom of the copper plate with 
У. To avoid shortcircuiting through the copper plate, 
the heater element and the bottom of the copper plate were 
painted with a high temperature nonconducting paint. То 
insure electric insulation, a 0.010 inch heat resistant 
gasket was placed between the heater and the copper plate. 

A one-quarter inch thick and seven inch square glass 
plate was used to cool the test liquid from the top. Five 
thermocouples with flattened beads were attached on the 
lower surface of the glass plate as shown in Figure 5. 

A one half inch wide plexiglass shim with a thickness 
of 0.125 inch was used as the spacer between the glass and 
the copper plate. The thermocouple leads from the glass 
plate were taken out through holes drilled into the plexi- 
glass shim as shown in Figure 6. 

The cooling chamber for the upper plate consisted of 
a half inch wide, one inch thick spacer with water inlet 
and outlet tubes on two sides as shown in Figure 7. Cooling 
water flowed through this space, bounded on bottom by the 
glass cooling plate and on top by a one-quarter inch thick 
and seven inch square plexiglass plate. To insure a 
watertight seal, two 0.0625 inch Neoprene rubber "O" rings 
were installed at each spacer-plate interface. 

To insure that the heat generated by the heater would 


go into the copper plate, a guard heater and two inches of 
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FIGURE 4. Heater Element 
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FIGURE 5. Location of the Thermocouples on the 
Glass Plate 
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FIGURE 6. Plexiglass Shim 
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FIGURE 7. Spacer of the Cooling Water Chamber 
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insulation was placed below the copper plate. The guard 
heater was constructed from nichrome wire in the same manner 
as the main heater element. Six thermocouples were installed 
in the insulation in two groups as shown in Figures 1 and 8. 
Power to the guard heater was adjusted so that the average 
temperature readings of the two groups were the same. Then 
the net heat loss through the insulation below the heating 
plate would be zero. This guard heater/insulation assembly 
was supported from the bottom with a one-quarter inch thick 
and seven inch square aluminum plate. 

The plates and the insulation were clamped together 
using four sections of aluminum "U" channel and four con- 
necting rods as shown in Figure 9. А leveling screw was 
attached to the bottom of each rod. By using the "U" channels 
to clamp the system together the load was distributed uni- 
formly. A one-quarter inch thick, half inch wide aluminum 
shim was placed on top of the plexiglass plate to distribute 
the loading to a larger area. A torque wrench was used to 
insure the same and desired amount of load was applied on 


each rod. 


CT INSTRUMENTATION 
1. Temperature Measurement 
The temperature distribution on the copper and glass 
plates was determined bv five ungrounded junction, copper- 
constantan thermocouples on each plate. These ten thermo- 


couples were wired into a thermocouple switching box, and 
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CURE 8. Location of Thermocouples in the Insulation 
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recorded separately by using a digital readout Numatron 
which displays the temperature reading in ?F directly. 

A thermocouple was attached to the inlet of the 
cooling water intake manifold and a second one to the outlet 
of the discharge manifold. To check these two temperature 
readings, an uncalibrated thermocouple was attached to every 
water inlet and outlet tube of the cooling chamber. These 
ten thermocouples were wired in parallel in two groups and 
each group was then wired into the switching box. These 
groups of five thermocouples in parallel gave average inlet 
and exit temperature readings of the cooling water. 

To control the heat loss through the insulation 
below the heater element, two groups of thermocouples were 
placed in the insulation. Each group contained three 
thermocouples and these six thermcouples were wired into 
the switching box separately. Аз mentioned before these 
thermocouples were used together with the guard heater to 
insure that the heat generated by the heater element was 
transferred into the test chamber through the copper plate. 

Thermocouple calibration procedure is given in 
Appendix A. 

2. Power to Heater Elements 

A Lambda regulated DC power supply was used to 
provide input voltage to the heater. А calibrated resistor 
was placed in series with the power supply and the heater. 
With this arrangement, voltage readings across the calibrated 


resistor and across the heater were taken. Knowing the 


59 





voltage of the calibrated resistance, the current through 

the circuit could be determined. The product of this current 
and the voltage across the heater gave the input power. By 
this arrangement the value of the input power could be 
accurately reset for different experiments. The voltage 
across the heater and across the calibrated resistance were 
determined using a Keithley 168 Autoranging DMM digital 
voltmeter. A sample calculation is provided in Appendix B. 

A variac was used to provide input voltage to the 
guard heater. The voltage across the guard heater was 
determined using the same digital voltmeter. 

TO ensure constant voltage supply, the input power 
for all electric equipment was taken from the output of an 
AC voltage regulator. 

3. Water Flow Rate into the Cooling Chamber 

A submersible electric pump was used to pump the 
cooling water from a 25 gallon reservoir. The flow rate 
was adjusted by controlling the voltage to the pump with a 
variac. A Fischer and Porter Co. flowmeter with a maximum 
rate of 0.6 GPM was used to determine the flow rate. The 
percentage of the maximum flow rate was marked on the 
flowmeter. 

A general arrangement of apparatus with the 


instrumentation is shown in the photograph in Figure 10. 
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ED. PROCEDURE 


A. APPARATUS ASSEMBLY 

The calibrated thermocouples were attached to the 
plates as shown in Figures 2 and 4. The plexiglass shim 
was placed on the glass plate and a watertight interface 
was accomplished by means of silicon rubber. A thick 
continuous silicon rubber film was laid on the upper sur- 
face of the plexiglass shim. The apparatus was then filled 
with the test liquid up to the edge of the silicon rubber 
layer. The next step was to attach the overflow tube which 
connected the test chamber to the expansion tank and to 
flll the whole tube with the test liquid. The copper plate 
was then placed over the glass plate and lightly pressed 
down. Once the excess liquid and silicon rubber flowed 
out, the test chamber was filled with the test liquid and 
all air bubbles were eliminated. Generally this result was 
accomplished by the third try. 

Assembly of the whole apparatus had to be accomplished 
quickly in order to avoid hardening of the silicon rubber 
film between the plates. Once the silicon rubber had 
hardened, a uniform thickness between plates could never 
be attained. 

саке from the bottom, first the aluminum support 
plate was placed on the two bottom "U" bars. The insulation 


with the guard heater and thermocouples was placed on the 
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support plate. The test chamber and the cooling water 
chamber were then placed on the insulation. The whole 
apparatus is shown in Figure 1. A torque wrench was used 
to compress the sections between the "U" bars. Starting 
with low torque values and increasing the load a foot- 
pound at a time an even load distribution and even silicon 
rubber thickness between the plates was insured. After 
measuring the distance between the copper plate and the 
glass plate at each corner with an inside micrometer, the 
test section was insulated on the sides. Later the apparatus 
was levelled using the leveling screws and a small bubble 


level on top of the apparatus. 


BEZZERROCEDURE 
l. Measurements 

Data of a usual run consisted of readings of the 
twentv thermocouples, voltages across the calibrated resistor 
and the heater elements and the cooling water flow rate. 
A complete experimental run usually lasted about three hours. 
A set of data was recorded every half hour. Steady state 
was considered established and final readings were made when 
the temperatures in the whole system varied less than 0.2 °F 
over half an hour. 

The recorded temperature data consisted of the 
inlet and exit temperatures of the cooling water, the temper- 
atures of the glass and the copper plates in the test chamber 


and the temperatures in the insulation below the test chamber. 
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Power to the heater element was monitored by the 
voltage across the heater. The voltage across the heater 
element was increased by two volts at a time starting at 
MITO TES, to a maximum of 25 volts. In the vicinity of the 
critical Rayleigh number the voltage was increased by smaller 
steps to obtain data at desired points. 

2. Determination of Dimensionless Numbers 

The data recorded at steady state included ten 
temperature readings for the copper and the glass plates 
of the test chamber. The average of the five temperature 
readings for each plate, gave the average temperature of 
the corresponding plate. The difference between the average 
temperatures of the two plates was defined as AT. А sample 
calculation is contained in Appendix B. 

All liquid properties were evaluated at the film 


temperature defined as 


The Rayleigh number was calculated using the distance 
between the two plates of the test section as the charac- 
teristic length. The Rayleigh number also can be written 


as the product of the Grashof number: 


g o AT 1? 


2 
V 


Gr = 
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ani the Prandtl number: 


"9 
= 
| 
o |< 


Voltages across the heater and the calibrated 
resistor were recorded for each run. Knowing the voltage 
across the calibrated resistance, the current through the 
heater circuit could be determined. The product of this 
current and the voltage across the heater gave the input 
power as shown in Appendix B. 

Heat leakage through the plexiglass shim between the 
two plates and heat leakage through the insulation below 
and on the sides of the test chamber were determined. The 
difference between the input power and the heat leakages 
gave the actual heat transferred to the test liquid which 
was defined as От and a sample calculation is contained in 
Appendix B. 

The Nusselt number was then calculated. By 
definition the Nusselt number is: 

ET, 
Nu < = 


Mino p 0 


Thermocouples located at the inlet and the exit 
manifolds of the cooling chamber gave the temperature rise 
of the cooling water in the cooling chamber. Together with 


the flow rate, this temperature difference, defined as 
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А gave a check on the calculated heat transfer into the 


Tew’ 
test chamber. A sample calculation is contained in Appendix 


в. 


ES TEST LIQUIDS 

A general idea of the required properties of the test 
liquid for the present design was obtained from the pre- 
liminary experiments conducted with distilled water. 

Finding the desired properties of the test liquids was 
one of the major difficulties encountered during this study. 
All desired properties of a selected liquid were not con- 
tained in one reference. In some cases the properties were 
given for a temperature range which did not include the 
temperatures encountered in the experiments. 

Three liquids were used.  Prandtl numbers varied from 
34 to 476 for the temperature range of the experiments. 
Properties of these liquids permitted temperature differ- 
ences in the range of 4 °F to 33 °F with very small layer 
depths. 

The first test liquid was a petroleum based oil commer- 
cailly known as Mobiltherm light 603. Properties were 
obtained from a technical bulletin [Ref. 14] published by 
the Mobil Research and Development Corporation. 

The second test liquid was ethylene glycol. Properties 
were obtained from the Handbook of Heat Transfer Media 


[Ref. 15] and Thermophysical Properties of Matter [Ref. 16]. 





The third test liquid was a glycerol-water solution. 
Composition of the solution was 56$ glycerol and 44$ 
distilled water by weight. Properties of the liquid were 


obtained from Glycerol [Ref. 17]. 
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IV. DISCUSSION AND CONCLUSIONS 


The hypothetical system of a fluid contained between 
two infinite, horizontal, conducting surfaces was closely 
approximated by containing the very thin liquid layer 
between two parallel square plates. Heat transfer through 
this very thin liquid layer bounded on top by a cooled 
glass plate and on bottom by a.heated copper plate was 
measured. Experiments were conducted for a range of Ray- 
leigh numbers from 350 to 4100 with three different liquids. 
The Prandtl numbers varied from 34 to 477 under these 
conditions. A summary of the data obtained from these 
experiments is presented in Tables I to IV. 

Table I contains the results of the experiments with 
Mobil 603. A plain glass plate was used as the cooled top 
surface. At high temperatures a short circuit between the 
heater and the copper plate was observed during these 
experiments. Later a very thin heat resistant gasket was 
placed between the heater and the copper plate. Data of 
these experiments is not used in obtaining the correlations. 

Table II contains results of experiments with Mobil 
603. During this second set of experiments a glass plate 
with its surface facing the copper plate coated with an 
electrically conducting transparent coating was used as 


the cooled top surface. 
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Tables III and IV present results of experiments with 
ethylene glycol and 56 wt $ glycerol-water solution respec- 
tively. A plain glass plate was used during both of these 
experiments as the cooled upper surface. 

The data under the LnRa and LnNu columns are plotted 
in Figure 13 as LnNu vs LnRa. The Nusselt number is about 
1 from a Rayleigh number of 350 to a range of Rayleigh 
numbers from 1500 to 1700, depending on the test liquid. 

Ra = 1600 + 100 was defined as the critical Rayleigh number. 
In their earlier experimental works Schmidt and Reiher 

[Ref. 9] and Ernst Schmidt and Silveston (Ref. 12] deter- 
mined critical Rayleigh numbers from 1600 to 1800 for 
pxous fluids. Up to this critical condition heat is 
transferred by conduction only. The deviation of Nusselt 
number in the conduction regime increases for lower tem- 
perature differences where the uncertainty is higher. 

At the Rayleigh numbers between 1500 and 1700 convection 
appears and the Nusselt number increases with the increasing 
Rayleigh number. The Nusselt number is nearly proportional 
to the Rayleigh numbers up to Rayleigh number of about 3000. 
A similar relationship was observed by Schmidt and Silveston 
[Ref. 12] for the same region which they called as the 
"Creeping convection" region. The authors obtained a 
correlation for the data in this region as: 


۲۲ 0 ОЕ 2 во 
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Correlations obtained for the data of this study are 
presented in Table V. 

A change in the slope above Ra = 3000 can be observed 
in Figure 13. The starting point of this change in the 
slope differed from one fluid to another. In the case 
of Glycerol-water solution it started well below Ra - 3000. 
Schmidt and Silveston [Ref. 12] concluded that this is the 
starting point of another mode of the convection heat 
transfer regime. 

The critical Rayleigh number which is the starting point 
of the convection heat transfer regime ranged from 1480 to 
1702. The critical Rayleigh number was determined by the 
intersection of the paired lines drawn through the data of 
conduction and convection regimes for each fluid in Figure 
13. These critical Rayleigh numbers and the corresponding 
Prandtl numbers are presented in Table VI and are plotted 
as LnRa-, vs LnPr in Figure 14. The critical Rayleigh 
number is apparently a function of the Prandtl number. 

The correlation obtained from this data is 
00 0 


Васр = 1103 (Pr) 
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V. RECOMMENDATIONS 


The uncertainty in the temperature difference was one 
of the biggest among the uncertainties. As the temperature 
difference increases the size of this uncertainty decreases. 
One could improve the precision of the temperature measure- 
ments by employing higher temperature differences. 

Another difficulty was attachment of the thermocouples 
to the plates of the test section. Even very small thickness 
of the flattened bead causes high uncertainty in the measure- 
ment of the distance between the two plates of the test 
chamber due to the very thin liquid layer thickness. 
Intrinsic thermocouples could be an answer to this problem. 

The corners of the plexiglass shim were chamfered 
to create a gap between the two plates of the test chamber. 
The distance between the two plates was determined by an 
inside micrometer from these gaps at four corners. In- 
creasing number of measurements could decrease the uncertainty 
of the Rayleigh number calculations. 

A flow visualization technique could be used to demon- 
strate visually the establishment of convection in a future 
send In tills area. Ernst Schmidt and Silveston [Ref. 12] 


used a shadowgraph technique and obtained good results. 
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APPENDIX A 


THERMOCOUPLE CALIBRATION PROCEDURE 


The accurate determination of the fluid and test section 
temperatures was an absolute necessity for this experiment. 
For this reason, precise calibration of the thermocouples 
was required. 

A Rosemount Calibration System, with a constant tem- 
perature oil bath, was used for the calibration. The 
eighteen thermocouples were suspended several inches into 
the oil bath. A Platinum Resistance Thermometer in conjunc- 
tion with a commutation bridge was used as a standard. The 
calibration was conducted over a range from 70 °F to 150 “Е. 
The maximum uncertainty for the standard thermometer for 
this temperature range was + 0.005 °F. 

The thermocouple readings during the experiments ranged 
from 74 °F to 130 °F. The maximum difference between the 
temperature reading of the standard thermometer and the 
thermocouples was 0.9 °F for this temperature range. For 
the same temperature range the maximum difference in tem- 
perature readings of the thermocouples was 0.1 °F. 

The calibration process was performed with the same 
switch-box and Numatron used during the experiments. In 
other words, there was no recording instrument changes 


once the system was calibrated. 
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A number was attached to every thermocouple after cali- 
bration. Positions of the thermocouples, except the ones 
used in the cooling water system, are shown in Figure 10. 
Thermocouples number 8 and 2 were located in the inlet and 
the exit cooling water manifolds respectively.  Thermo- 
couples attached to the inlet and the exit tubes of the 


cooling water spacer were numbered 19 and 20 respectively. 
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APPENDIX B 
SAMPLE CALCULATIONS 


TO represent a sample calculation, the data recorded 
on 23 July 1977 is used. The test liquid was glycerol- 
water solution. Composition of the solution was 56$ 
glycerol and 44$ distilled water by weight. 

Steady state was reached in two hours twenty minutes. 
Maximum temperature change in fifty minutes was less than 
0.2 °F for the copper and the glass plates of the test 
chamber. The temperature change for the cooling water 
was less than 0.1 °F for the same period of time. 

A sketch of the control volume for the energy balance 
on the test chamber indicating the major heat transfer 
components involved is shown in Figure 12. 

Sample calculations of the supplied power (Op) y the 
heat leakage through the plexiglass shim (Qr) , the heat 


loss through the insulation (0 ), the heat transferred to 


INS 


the cooling chamber (Q...), the Rayleigh number (Ra), the 


CW 
Nusselt number (Nu) and the Prandtl number (Pr) are given 
below. 

The heat loss through the insulation consisted of the 
heat loss through the insulation below the test chamber 
(Од) and the heat loss through the insulation on the sides 


of the test chamber (Q.). The sum of the heat losses from 


5 
the sides for each plate of the test chamber was equal to the 
heat loss for an equivalent cross sectional area at the film 


temperature. 
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FIGURE 12. Energy Balance in Test Chamber 





SAMPLE CALCULATIONS 


A. DATA 
ШЕ Т = 93,9 °F 
“ы 2 = 73.1 ФЕ 
Т аз 
Ш. 4 = SG TF 
М5 - 58177 
UNO = 88.3 °F 
H7 = 87.9 °F 
Е В = 72.8 °F 
9 = 80.6 °F 
ПЕ ЩО = 81.1 °F 
peu) Е СОЕ 
Ш 12 - 81.3 °F 
ШЫ Зз = 81.37 F 
Mee 14 = 94,1 °F 
سا‎ 94 Ser 
Дар 16 = 9998 °F 
Цас іс = 04% er 
0-7 
oe | 7 
eco = زور‎ Е 
Е - 13.98 V 
En = 2.31 V 
E = 4.10 V 





R = OLS 
CW 2 120.2 lbm/Hr 

p = 75 ФЕ 

Кр - ПИОИВЕЦ НЕ /ЕЕ/29Ғ 

А = 0.0903 #+2 

P 

I = 0.01167 ЕЕ 

Kons = O 096 I BtuZHr/£ft/°F 
P 2 

Ap = оо ‘ЕЕ 

L5 = 0.02333 ft 

Lo = 0.02983 ft 

Ceu = 0.998 Btu/Lbm/°F 

8 ےا تد‎ 0 ^ 1/*r 

V = ZO x l £t^/sec 

Kr ro = O22 Btu Hr/fE SF 
= 3 

PLIO Е Ом ЕЕ 

CLIQ = 0.7636 Btu/Lbm/°F 
2 2 

Ar ro = 0.25 ft 


TEMPERATURE CALCULATIONS 


E eee melo oe tear ote. Temperature (Тс) 


Тс 3 + Тс 4 + Тс 5 + Тс 6 + Тс 7 


C 5 


OT 555 LTE 98.3 + 87.9 
5 


Ber 0 22“ 
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AueragenGiass Plate Temperature (Тс) 


Тс 9 + Тс 10 + Тс 11 + Тс 12 + Тс 13 


С 5 


80.6 + 81.1 + 80.9 + 81.3 + 81.3 
Б 


Saucen E 


Temperature Difference Between the Two Plates (AT) 


дар 


II 
H 
۱ 
на 


= == = ° 
C G С Be ح7۰۳7‎ OG 


) 


Film Temperature (To 


) 


Temperature Rise in the Cooling Water (AT Arg 


= - - = = £ 
AT ow Tow 2 Тси 8 psc 275 Е 


еее рое Dil rerence іп the Insulation (AT ug) 


" И Тус J + Tes Л Tos 18 В Тус 14 + Тус hee Туус 


INS 3 В 


INR 93.7 9421 + 94.3 + 93.8 
7 3 a کچ‎ 


SO = 091 0.11 £ 0:3 °F 


oul 





7. The Difference Between the Room Temperature and 
tee tin Temperature (АТ ) 


peed 84.5 - 75.0 = 9.5 °F 


C. POWER CALCULATIONS 


1. Supplied Power (Q5) 


Е یی‎ 


Ор = В 2.031 


15:90 WN = 54.3 BtEn/hr 


2. Heat Leakage Through the Plexiglass Shim (Q,) 


Ве © 120) 00903) (7:1) 


ІС ут 0.01167 


6.6 Btu/hr 


3. Heat Loss Through the Insulation (Qing) 


Qins = 9p * Ys 


ENSE 7 | pice cles 


Рв HS 


. (0.096) (0.25) (0.1) , (0.096) (0.1244) (9.5) 
0.0833 0.02083 


z20019)901:105551— 75.53 Btu/hr. 
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4. Heat Transferred to the Test Chamber (Qn) 


Т 


Btu/hr‏ 42.2 = ر6 5ت 5۔5۰۱ 


5. Heat Transferred to the Cooling Chamber ) 


Cow 


Jew ^ Few “си Stew 


(ШОО 2y = 42.0 Btu/hr 


|| 


ши SIMENSTONLESS NUMBERS 


аео Number (Ва) 


быр Т S 


a V 


Ra = 


4 23 


٠٠1-03 1:167 0٦ 
کر ...بت‎ 


(29174) (2.683 x 10 
00032 | 
(71.02) (0.7636) (3600) 


192032 
2. Nusselt Number (Nu) 


Nu = 


ОО ох) 


10.232) 17.02) (0.25) = 1.2106 
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3. Prandtl Number (Pr) 


x 
Pr = — 
а 


o gea 
АСО 7656) 


47.6507 
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APPENDIX C 


UNCERTAINTY ANALYSIS 


The uncertainties for the variables and the dimension- 
less numbers in this experimental study were calculated 
by the method proposed by Kline and McClintock [Ref. 18]. 

The second-power equation of Kline and McClintock was 
used for the calculation of uncertainties in the values 
obtained experimentally. 

The basis for uncertainties in the temperature readings 
was the calibration of the thermocouples. For all measured 
quantities the accuracy of the measuring instrument was the 
basis for uncertainties. Uncertainties for the properties 
of the test liquids were given in the references where these 
properties were obtained. 

As an example, the calculation of the uncertainty in 
the Rayleigh number for the same case taken for sample 
calculations is given below. The Rayleigh number was 


defined as: 


g B ATIL 
va 


Ra = 


and the uncertainty was calculated by the seoncd-power 


equation of Kline and McClintock as: 
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dRa _ 
Ra 


ав 


13558 


2 cim 2 


)^ € (XD € 302 € (3) + ) 


ДИ 


(0.0121)2 + (0.085) 


dL, 2 du 2 , (da, 2 


V а 


2 2 


+ 3(0.060)% + (0.001)% + (0.0102)? 


Using the same formula uncertainty in the Nusselt number 


and the Prandtl number was found to be 16.4% and 1.143% 


respectively. 


are 


the 


Me Values of the uncertainties for other variables 


listed below. Because of 


the low temperature difference, 


uncertainties are higher than average for this case. 


Quantity 


dor to, Prio 


d А 


8 LIO 


4 rg CLIO 


d, IQ LIO 


ак,/%ь 


ак k 


Ins’ INS 


66 


Uncertainty 


0 0010 


0.0100 


ШЕШІ? 


972 00/06 


02.058 


8,۷۸" 0 








Quantity 


dL/L 


dL A. 


INS' INS 


dE, /E,, 
dE,/E, 
dR/R 
dAT/AT 
d8/8 
۹٥+ 

30, /0, 
dOrus/ 1s 


20/0. 


da/a 


67 


Uncertainty 


0.060 


0.050 


0.0010 


0.0050 


0.0005 


0.0850 


000 


0.0051 


021201 


0. 7 


0. 1255 


0.0102 





Quantity Uncertainty 


ЧАТ то/Атто 0.020 
dAL/Az 0.020 
ЧА مم‎ 0.020 
дА А, 0.020 
am Ly 0.050 
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